function DecDev=SubFun_UnitIRFDec(PP,SS,MODEL,VarList,InputDev)
%% Baseline Input Deviation (all 0s)
FullVarList     =   fieldnames(InputDev);
InputDev_0      =   struct();
for ii=1:length(FullVarList)
    TempVar         =   FullVarList{ii};
    InputDev_0.(TempVar) ...
                    =   zeros(size(InputDev.(TempVar)));
end

%% Individual Effects
OutputDev       =   struct();
for ii=1:length(VarList)
    TempVar         =   VarList{ii};
    TempInputDev    =   InputDev_0;
    TempInputDev.(TempVar) ...
                    =   InputDev.(TempVar);
    OutputDev.(TempVar) ...
                    =   SubFun_Solver_CPaths(PP,SS,MODEL,TempInputDev);

end

DecDev      =   OutputDev;
